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ABSTRACT 

We extend the work by Appl and Camenzind (1988) for special relativistic magnetohydro- 
dynamic (MHD) jets, to fully general relativistic studies of the standing shock formation for 
accreting MHD plasma in a rotating, stationary and axisymmetric black hole magnetosphere. 
All the postshock physical quantities are expressed in terms of the relativistic compression ratio, 
which can be obtained in terms of preshock quantities. Then, the downstream state of a shocked 
plasma is determined by the upstream state of the accreting plasma. In this paper sample so- 
lutions are presented for slow magnetosonic shocks for accreting flows in the equatorial plane. 
We find that some properties of the slow magnetosonic shock for the rotating magnetosphere 
can behave like a fast magnetosonic shock. In fact, it is confirmed that in the limit of weak 
gravity for the upstream non-rotating accretion plasma where the magnetic field lines are leading 
and rotating, our results are very similar to the fast magnetosonic shock solution by Appl and 
Camenzind (1988). However, we find that the situation becomes far more complicated due to the 
effects of strong gravity and rotation, such as the frame dragging-effects. We show the tendency 
that the large spin of the black hole makes the slow magnetosonic shock strong for the accretion 
solutions with the same energy-flux. 

Subject headings: accretion — black hole physics — MHD — relativity — shock waves 



1. Introduction 

Recent theoretical and observational developments suggest that various members of active galactic nu- 
clei (AGN) involve supermassive black holes (e.g., Rees (1997)). According to the 'standard unified scenario' 
an AGN reveals itself as a quasar or a Seyfert nucleus during the earlier stages when the accretion rate is 
high and the black hole is spinning up, while in later stages when the accretion slows down and the black hole 
starts losing its rotational energy by spinning down it evolves to a radio galaxy. In the past years physics of 
black hole magnetospheres has been extensively studied, mainly in connection with winds, jets, and energy 
extraction from black holes in rotation-powered AGN (Blandford & Znajek (1977); Znajek (1977); Thorne, 
Price, & Macdonald (1986); see, e.g., Begelman, Blandford & Rees (1984) for comprehensive review). An im- 
portant step for studying the magnetosphere in the Kerr background space-time was originated by Blandford 
& Znajek (1977) in the force-free limit. The next logical step will be to extend these magnetospheric studies 
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to accretion-powered AGN — especially to Seyfert nuclei where the accretion rate should be relatively mod- 
erate. Phinney (1983) first extended his general studies of black hole magnetospheres to accretion-powered 
AGN, adopting the magnetohydrodynamics (MHD). In the early 1990s, Takahashi et al. (1990, hereafter 
TNTT90), Nitta, Takahashi & Tomimatsu (1991) and Hirotani et al. (1992) investigated both inflow and 
outflow of material in black hole magnetospheres. Their approach was a general relativistic unified MHD 
description of both matter accretion process and electromagnetic process in the magnctosphere. Relaxing 
the force- free limit, these authors explored the roles of interactions between the accreting fluid matter and 
electromagnetic fields. 

Recently, ASCA observations indicated evidence that the Fe lines observed from some Seyfert nuclei 
are emitted from regions very close to the central black hole (e.g., Tanaka et al. (1995); Nandra et al. 
(1997)). Iwasawa et al. (1996) interpreted the behavior of Seyfert MCG 6-30-15 deduced from the long 
ASCA observation as an indication that the black hole is rotating extremely fast, near the maximum limit. 
Reynolds & Begelman (1997) pointed out that such extremely fast rotation contradicts with the standard 
unified scenario. Instead, these authors showed that if an X-ray point source somewhere above the hole on 
the rotation axis illuminates the infalling gas within the inner accretion disk radius, a Schwarzschild hole 
is consistent with the observation. However, it is not clear how such an X-ray source can be created in the 
specified location. This issue clearly points to the importance of investigating, for accretion-powered AGN 
also, the basic physics of the vicinity very close to a black hole — especially the regions between the inner 
boundary of the accretion disk and the event horizon. As our first step toward such investigation, therefore, 
we explored the standing shock formation in accreting plasmas in black hole magnetospheres. These studies 
may shed valuable insight to problems such as how X-ray sources can be created near the event horizon. 

We consider the central engine of an AGN as a black hole magnetosphere, where magnetized plasma 
surrounds a black hole and infalling accreting flows and outgoing wind /jets would be generated from the 
surrounding plasma. We formulate the MHD shock conditions in Kerr geometry for such plasma. The 
formation of shocks is based on the existence of multi-magnetosonic points in the accretion solution, while 
for the outgoing plasma the sub-magnetosonic solution is allowable at distant regions. This is because, for 
instance, the accreting flows initially ejected from a plasma source with low velocity must be terminally 
super-fast magnetosonic at the event horizon. At the shock front, the flow transits from super-magnetosonic 
to sub-magnetosonic, so that the accreting flows with a shock must pass through a magnetosonic point on 
each side of the shock front. This situation is quite similar to the case of hydrodynamical accretion onto a 
black hole (e.g., Chakrabarti (1990a); Sponholz & Molteni (1994); Lu et al. (1997); Lu & Yuan (1998)). The 
trans-magnetosonic MHD flow solution was discussed by Takahashi (2000a, b, 2001). Along the magnetic 
field line, the five physical quantifies are conserved: the total energy E and angular momentum L, the 
angular frequency of the magnetic field line £If, the particle number flux per magnetic flux tube ry, and 
entropy S (see, e.g., Camenzind (1986a)). When these conserved quantities are specified at the plasma 
source, the location of the fast/slow magnetosonic points and the Alfven point are determined. Takahashi 
(2000a) obtained multi-magnetosonic point solutions and found two regimes of accretion flows - "hydro- 
like" and "magneto- like" . The hydro-like accretion would transit to magneto-like accretion by the shock 
formation. However, we will postpone, until later, the problem of joining these two types of solutions by 
shock formation. The reason is that in order to do so we will have to carry out detailed parameter search for 
trans-magnetosonic MHD flows, but then we will have five field-aligned physical quantities in both upstream 
and downstream solutions, not a trivial situation. This problem will be additional in a subsequent paper 
(Rilctt et al., in preparation). This paper is meant as a starting point for our long-range investigation of 
shock conditions for accretion flows in the Kerr geometry. Therefore, here we will only solve the cold trans- 
fast MHD equations for upstream accretion and discuss the shock properties at the shock fronts. At the 
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shock front the plasma would be heated up. So, the postshock accretion should be treated by a hot MHD 
accretion model (Takahashi 2001). In this paper we will not yet solve explicitly the hot trans-magnetosonic 
solutions for postshock accretion but instead treat the shock front location as a free parameter. However, by 
joining preshock and postshock solutions, the shock location would be determined at one (or more) place(s). 

The main purpose of this paper is to explore the effects of rotation and general relativity on the MHD 
shock conditions for accreting plasma in a black hole magnetosphere. In order to do so, first in Section 
2 we present the basic equations for MHD accretion in the Kerr geometry. Then, in Section 3 we extend 
the work for special relativistic MHD jets by Appl & Camenzind (1988, hereafter AC88), by deriving the 
shock conditions for general relativistic MHD accretion onto a Kerr black hole. (Non-relativistic MHD 
outflow solutions with standing shocks were discussed by Chakrabarti (1990b).) Following AC88, all of our 
postshock physical quantities are expressed in terms of the relativistic compression ratio £. This compression 
ratio is the solution of a polynomial of eighth degree. Due to the additional factors, namely magnetosphere 
rotation and general relativistic effects, the mechanism of solution is far more complicated and tedious 
than in AC88. Our equations reduce to the AC88 equations in the weak gravity limit without plasma 
rotation. TNTT90 extensively studied MHD accretion flows in the Kerr geometry. As the next step our 
shock conditions are applied to the accreting MHD plasma flows as described by TNTT90. In Section 4 we 
present some examples of representative physically relevant shock solutions found for acceptable accretion 
flows onto the event horizon. Our results are presented for the equatorial flows and slow magnetosonic shocks. 
The "switch-off" shock (switching-off of the magnetic field at the shock front) and comparison between AC88 
and our results are also discussed. More thorough presentation of the results and extensive discussion of the 
physical implication of these results will be given in Paper II (Rilett et al., in preparation). The general 
cases including the fast magnetosonic shocks will be presented in Paper III (Rilett et al., in preparation). 
The possible application of our current work to some astrophysical problems, such as the Seyfert Fe lines, 
will be given in subsequent papers. Summary and conclusion are given in the last section 5. 



2. Basic Equations of General Relativistic Plasma Flow 

In this section we summarize the general relativistic MHD flows (see, e.g., Camenzind (1986b, 1989); 
TNTT90). We assume a stationary and axisymmetric magnetosphere and ignore its self-gravity. We also 
require infinite conductivity for the plasma flow. The background metric is given by the Boyer-Lindquist 
coordinates with the c = G = 1 units 

ds > = (i-gg:)^ + 4amr s sin2 % w a) 

Zj Za 

where A = r 2 — 2mr + a 2 , £ = r 2 + a 2 cos 2 6, A = (r 2 + a 2 ) 2 — a 2 A sin 2 9 and m and a denote the mass and 
angular momentum per unit mass of the black hole, respectively. The basic equations of relativistic plasma 
flows are as follows: (i) the particle number conservation 

K);a - , (2) 



where n is the proper particle number density and u a is the fluid 4-velocity; (ii) the conservation of total 
energy and momentum 



T a f = , (3) 
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where the energy-momentum tensor is given by 

T af3 = n^u a u f3 - Pg af3 + ^- (f^F^ + \g aP F 2 ^j , (4) 

fi is the relativistic specific enthalpy, P is the pressure, F a/3 is the electromagnetic field tensor and F 2 = 
F a PF al3 ; and (iii) the MHD condition 

u fi F aP = . (5) 
The magnetic field and electric field seen by a distant observer are defined as 

B a = \^M^ iF1& , (6) 
E a = F aP kP , (7) 

where k a — (1,0,0,0) is the time-like killing vector and r? Q( 3 7 5 = \f—g ^ap-yS- Here, we can introduce the 
angular velocity of the magnetosphere (see, Bekenstein & Oron (1978)) 

^--# L = -#^' (8) 

which is conserved along magnetic field lines. Then, we obtain the following useful expressions 

Ee = -V=dn F B r /G t , (9) 
E r = ^n F B e /G t , (10) 

where G t = g u + Qt^F- 

From the conservation of total energy and momentum (3) and the Killing equation + Xv^ = 0, 
where 1S a Killing vector, it follows that {Xfj,T fll/ )- l , — (see, Blandford & Znajek (1977)). Thus, for any 
stationary and axisymmetric system, we define the conserved energy flux 

= T^K = T? = (T t ") em + (T/Vid (11) 

and angular momentum flux 

= T^m v = T; = (T£) em + (T£) fluid , (12) 

where m v is the axial Killing vector with Boycr-Lindquist component (0,0,0,1), and the electromagnetic part 
and fluid part are labeled by "em" and "fluid" , respectively. The fluid and electromagnetic parts of the 
radial component of energy flux £ r are 

£fluid = (TDfluid = nfiu r u t , (13) 

and the fluid and electromagnetic parts of the radial component of angular momentum flux C r are 

-£fluid = ( r 0)fi"id = n^u r U(t> , (15) 

~ £om = (T ^ )cm ~I^l7^ + J ~^g7' (16) 
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where we have used the relation (9). 

We can express the radial energy and angular momentum fluxes as £ r = nu r E and C r — nu r L, 
respectively, where E and L are the total energy and the total angular momentum of the MHD flow seen by 
a distant observer and defined by 

E s (IT, 

L S itr,' <18) 

and 77 is the particle number flux per magnetic flux tube defined by 

n=-%rG t . (19) 

These quantities E, L and 77 are conserved along stream lines, which coincide with magnetic field lines 
(see, Camenzind (1986a)). 

The energy and angular momentum of the fluid are 

M 2 E - eG t 

M«t = M2 _ a > (20) 
M 2 L + eG^ 

W = - M"-a ' (21) 
where M 2 = 4nr] 2 (fj, / n) is the Alfven Mach number, e = E — QpL, = g t( p + g^^F and a = G t + G^Qf- 
The poloidal equation (the relativistic Bernoulli equation) of the magnetized flow is written as 

(l + u 2 p ) = (E/tf[(a-2M 2 )f-k} , (22) 

where 

k = -{g tt - 2g t4 'L + g^L 2 ) , (24) 

Up is the poloidal velocity defined by u 2 = —u a ua {A = r, 9), f? w = g 2 ^ — gttg<p<p — A sin 2 6 and L = L/E. 
We set Up > (u r < 0) for ingoing flows. When we neglect dissipative effects and cooling in the plasma, the 
flow is adiabatic with a specific enthalpy of the form (see, Camenzind (1987)) 

r p 

a* = m p + > ( 25 ) 

where m p is the rest mass of the particle. The toroidal component of the magnetic field can be expressed as 

G^E + GtL 

B * = ^ M 2 -a ■ (26) 

Figure 1 illustrates solutions for the poloidal flow equation (22) for a radial field geometry in the 
equatorial plane, where the plasma is cold (P = 0) with given parameters E, L and flp in the Schwarzschild 
spacetime. [Note that to determine the magnetic field configuration we should treat the equation of cross- 
field momentum balance (Grad-Shafranov equation); the general relativistic extension was obtained by Nitta, 
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Takahashi & Tomimatsu (1991). However, the task is very complicated, so that here we only use a given 
magnetic field configuration (see section 5).] The flow must pass through the fast magnetosonic point 
(TNTT90). That is possible only if there is some specific relation r\ = r](E, L,Qp); Only for the heavy solid 
curve with the arrow is this relation satisfied. The physical accretion solution starts from r = r; n j with zero 
poloidal velocity and reaches the event horizon r = rjj with a non-zero poloidal velocity. A slow magnetosonic 
shock should be located between the plasma injection radius r = n n j and the Alfven radius r — ta, if the 
shock conditions are satisfied. After the slow magnetosonic shock, the heated postshock flow falls into the 
black hole after passing the slow magnetosonic point, the Alfven point and the fast magnetosonic point, in 
this order. The radius of the light surface r — tl and the Alfven radius r — r\ do not change after the shock 
formation (see the next section). 



In order to explore the properties of MHD shocks associated with accretion flows near a black hole, we 
derive in this section the general relativistic version of MHD shock conditions — an extension of the work by 
AC88. Figure 2 shows a schematic picture of general accretion inflow from a plasma source, through a shock 
front, and onto a black hole. In this picture, the accretion originates from the plasma source (which can be 
the surface of a torus) rotating around a black hole. Strong shocks may be produced somewhere between 
the plasma source and the event horizon. 



In a complete solution of the accretion that includes a shock, the flow must satisfy a set of conditions 
on either side of the discontinuity. The jump conditions for arbitrary shocks in a relativistic MHD flow are 



where n a = (no,n) and the square brackets denote the difference between the values of a quantity on the 
two sides of the shock. 

We assume that the downstream flow velocity is radial (normal to the event horizon) and the shock 
front is perpendicular to the downstream flow n=(l,0,0) as shown in Figure 2. Then, we set 



3. MHD Shocks in Kerr geometry 



3.1. The Jump Conditions 



(e.g., AC88): 



[nu a ]n a = , 
[T afs ]n a = , 
[E] x n = , 
[B] • n = , 



the particle number conservation 

the energy momentum conservation 

the continuity relations for the electric field 

the continuity relations for the magnetic field 



(27) 
(28) 
(29) 
(30) 




{u\,u r 2 ,Q,ut) , 
(B r , B{, Bf) , 

(E{,E e ,Q) , 



(31) 
(32) 
(33) 
(34) 
(35) 
(36) 
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where equations (29) and (30) have been used. The subscripts "1" and "2" denote the preshock and the 
postshock quantities, respectively. From equation (9), we see that fi_F does not change across the shock. 
Equations (27) and (28) evaluated in the shock rest frame yield the following relations 

uyu\ — n 2 u 2 , (37) 

mmivriir)! -Pi + -^—{-E[E rl + B\B B1 + BfB^) 
oirgtt 

= n 2 ^ 2 {u r u r ) 2 -P 2 + -^—(BtB^), (38) 
oirgtt 

ni/iiuju? - — !— (B r B? + E{E e ) = , (39) 

r><fi r>r r>4> or 

nxli\u\u\ - = n 2 fi 2 u 2 u% - , (40) 

47r# tt ^gtt 

n^ulu\ - - J- ( + B[bA = n 2 ^ 2 u\ -J— ( + . (41) 



From the MHD conditions, we also obtain 



Eriu tl + J—g(Btu\ - B\u\) = , (42) 
{B r u f _ B f u;) J_ = _J| = (iru * _ B 2 ^) J_ , (43) 

u ti \J-g u t2 

(g^Ee + J~gB r )u{ + (g tl pE rl - ^~gB\)u\ = , (44) 
E rl u[ + E ei u{ = . (45) 

From the shock condition of (T r *)i = (T r *) 2 , we obtain 

(E - wi)! = (E - uL) 2 , (46) 

where ui = —gt<)>/ 9<t><j> is the angular velocity of the zero angular momentum observer (ZAMO) with respect 
to a distant observer, and from the shock condition of (T r *)i = (T r<t> ) 2 , we also obtain 

(g t4> E + g tt L) 1 = (g t4> E + g tt L) 2 . (47) 

From these relations, we see that E\ — E 2 and L\ = L 2 . Furthermore, from the particle number conservation 
and the continuity for electric and magnetic fields at the shock front, we also obtain the relations (fip)i = 
(CIf)2 and r/i = f] 2 . Thus, the field-aligned flow parameters E, L, tip and ij are conserved across the shock 
front. This means that the location of the light surface and the Alfven radius do not change by a shock 
generation. The locations of light surfaces r = n, are given by a — and the Alfven radius r = r\ is given 
by L — —{G^/G-t)\. Because of the entropy generation (S 2 > Si), the fast/slow magnetosonic points appear 
at different locations in the preshock and postshock flow solutions. The increasing entropy makes the Alfven 
wave speed at the Alfven point ma [= maw(^a)] in a postshock solution smaller than that in a preshock 
solution. 
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3.2. Dimensionless Parameters and Their Relations 

From the energy momentum conservation at the shock front, we have (T rt )i — (T rt )2, where T rt can 
be reduced to 



T rt = ff"(4 r uid + ^ m )-5 t0 (Auid + ^e m ) (48) 

(49) 



9* 



r B^B r (Q F -cj) 
Tijiu u t (l — ujI) + ■ 



and I = —Utjt/ut is the specific angular momentum of the plasma. Here, we define dimensionless parameters. 
First, we define the "magnetization parameter", which denotes the ratio of the Poynting flux to the total 
mass-energy flux seen by ZAMO, 

(t r - bjL r ) fluid 

then we can express equation (48) as 

T rt = g tt (£ r - ^)fl uid (l + a) = (nu r )(0u')(l + a) . (51) 
Then, the magnetization parameter can be reduced to 



Airri^u 1 ) pi 

By using (20) and (21), a is denoted in terms of M 2 , E, L and Qp, as 



(52) 



(GtpE + G , t £)G e~ha 
P 2 w e + {guE + gt4,L)M 2 ~ e - hM 2 ' ( ' 



where h = g tt (E — uiL). We can also express M 2 in terms of a, E, L and tip. 
Next, we define the following dimensionless parameters 



B% M- 



" S if = W~a • (54 » 

where q is the amplification factor for the toroidal magnetic field and £ is the shock frame compression ratio. 
From equations (37), (41), (43) we obtain the following relation 

1 - C = o-i (« - 1) • (57) 
The normalization of the 4-velocity gives an equation for u\ 



u\ = 



2 ( u % 1 

9tt + 2gt4>M 2 + g<p<pii 2 + 
where = vfi ' ju l is the angular velocity of the fluid, and 



(58) 



O _ eO f - p w 2 {gt4>E + gttL)M 2 

" 2 T^hMl ' (59) 
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When other postshock quantities £ or Mf are determined, u\ and Vt 2 can be obtained. 
From equation (38), we obtain 

1 _ _ _ + (--E r ^r + + B^B^ {B+B+h = Q 

£ ni,ui(u r u r )i 87T5 tt ni^i(u r u r )i 



£ r £ r + B B g = ^—B B e , (61) 
G f 



Using the relations 

3 u _ 9tt® D g 
I 2 

„2 



we get 



with 



)i-(B S ) 2 = (B^)i(l-O. (62) 



1 - | = n - ^ + (rj 2 - l)7i (63) 



n = Pl ,f\ = l-j + ^-(q-l) + Xi -(q 2 -m , (64) 



* = 2B^Mp (65) 



Tl = tU-j^T' (66) 

which depend on £, g and upstream parameters M\ ( or o\ ), £7, L and ftp. The relation between M| and 
Mf are given by the poloidal equation (22). (At the shock location r = r sh , both M 2 = Mf and M 2 = Mf 
are solutions of the poloidal equation.) 

In the following, for preshock accretion flows we restrict ourselves to the cold limit (Pi = 0). Using the 
definition of £ and the equation of state for a Boltzmann gas with the polytropic index T 



Mi = m p (67) 

r p 2 " 

r - 1 n 2 m P/ 



= «p(1 + A^) - (68) 



we find with equation (64) 



Combining equations (54), (56) and (57) gives the quadratic equation for £ 

c 2 - 



. Gt ciClp 



C+(-^-h^Vk = o, 



Mf n F - 



where we use the relation 51 — Hp = aM 2 /G<p. We are now able to eliminate u\ and q from equation (69). 
After considerable manipulations we obtained a polynomial of eighth degree in £ 

8 

Y,c i (M 2 ,r:m,a,E,L,n F ,T 1 )e = . (71) 

i=0 
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The coefficients a are dependent only on upstream parameters, except for T which is a function of the 
downstream temperature involving the modified Bessel functions (Synge 1957) 



r(e) = 1 



1 fK^l/Q) 



e \k 2 (i/o) 



1+3 



with 6 = kT/m,p. 6 is related to £ though 

9 = -grrUluiiul/ulfll/Z . 



(72) 



(73) 



Thus, the compression ratio £ is the solution of a polynomial (71), which has to be solved simultaneously with 
equation (72) for the polytropic index T for the shocked plasma. The concrete expressions of the coefficients 
Ci are quite lengthly, so that we only show numerical results in the next section. This polynomial of eighth 
degree has in general several real solutions corresponding to the different shock transitions. The downstream 



quantities (, q, 
respectively. 



£1 2 and II are also obtained from equations (70), (54), (58), (37), (59) and (64), 



4. Results and Discussions for Slow Magnetosonic Shocks 

To see the general shock behavior which depends on flow parameters with the hole's spin a, we present 
solutions to the coupled equations (71) and (72) for a range of black hole spins. For computational reasons, 
instead of equation (72) we use a polynomial approximation given by Service (1986). In this paper, we 
restrict ourselves to slow magnetosonic shocks, which should be located between the plasma source and 
the Alfven point. Then, we will only treat the sub-Alfvenic region of M < Ma on a trans-fast MHD 
accretion solution, which is considered as a preshock accretion solution and is assumed cold. This preshock 
accretion is super-slow magnetosonic. The postshock accretion would be heated up at the shock front, and 
the slow magnetosonic wave in the hot plasma has a non-zero wave speed (usm)2 > 0, where usm is the 
slow magnetosonic wave speed (see, TNTT90). Then, the hot postshock accretion must satisfy the condition 
< u P 2 < («sm)2 at the shock location; that is, the postshock flow must be sub-slow magnetosonic. All 
our preshock accretion solutions presented in this paper refer to physically acceptable accretion flows, in the 
sense that the critical condition (as shown as an example in Fig. 1) is satisfied. 

Since we have four conserved parameters, E, L, Qp and 77, two of these quantities are kept the same 
for all spin parameters. (The physical descriptions of the parameters are summarized in Table 1.) Fixing 
L and Qp would prevent the formation of shocks across the entire spin range, and so E and r\ were held 
constant. Then, physical parameter sets for the trans-fast MHD accretion solutions are calculated. The 
trans-fast MHD accretion solutions were then obtained by letting VLp increase and L (and thus e) decrease 
linearly with increasing spin a (Table 2). For simplicity, we set the flow to be in the equatorial plane of the 
radial field geometry B r = i?i n jG t /S and B 6 = (and then u e = 0). B; n j is the strength of the field at the 
plasma injection point. Here, we should note that the radii of r = rn = m + \/m 2 — a 2 , r = rp, r = tl, 
r = ta and r = ri n j arc different for each flow solutions. For larger spin a, the region of ta < r < ri n j, where 
the slow magnetosonic shock is expected, becomes narrow and shifts inward (toward the horizon). 

In figures 3 and 4 we show the effect of black hole rotation on the nature of shock and accretion flows. 
The spin parameters selected are a/m = 0.95, 0.3, 0.0, —0.2 and —0.4. (In this paper we only calculate the 
cases SIf > Wh-) In the figures the solid curves are for the maximum a (= 0.95m), and the spin decreases as 
it changes from the solid to the dot-dashed curves which are the minimum at a = —0.4m. Positive a refers 
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to black hole and magnctosphcrc corotation, while negative a refers to the counter-rotation. No solutions 
were found for a < —0.4m. This is due to the choice of conserved quantities. 

Figures 3a, 3b, 3c and 3d show the shock frame compression ratio £, thermal pressure II, polytropic 
index T, and temperature parameter 9, respectively, as a function of v\, where v r = — {-/A/ /S)u r /u l is the 
positive radial 3-velocity of the infalling plasma as seen by a ZAMO. The parameters II, Y and 6 refer to 
the downstream values. Figure 3a shows in general that the compression ratio £ increases and then decreases 
with velocity. Figures 3a and 3b show that both £ and II increase with increasing black hole rotation. 
Thermal energy kT (or 9) and polytropic index Y for postshock flow have a weak dependence on the black 
hole spin, but the general trend is that for a given velocity, Y slightly decreases while 9 slightly increases with 
increasing \a\. Figures 4a, 4b and 4c show the effect of velocity on three other interesting parameters — the 
ratio of downstream toroidal magnetic field to upstream toroidal magnetic field q, the ratio of downstream 
number density to upstream number density n^/ni, and the upstream magnetization parameter u\. We find 
that q generally decreases while the ratio of the postshock to preshock density increases with the velocity. 
The effect of larger spin is to enhance the density jump while the azimuthal magnetic field amplification and 
magnetization are suppressed. In Figure 4a, we can see the switch-off shock, except for the counter-rotating 
black hole cases. The switch-off shock occurs in the limit when the flow velocity equals the Alfven wave 
speed at r s h = ta- At this limit, for the postshock flow the toroidal magnetic field Bf seen by distant 
observers becomes zero (q = 0). When the shock happens close to the Alfven point, a very hot plasma region 
is generated (the small spots in Fig. 3d). This region could emit X-rays and 7-rays around the black hole. 

In the Newtonian analogy, we expect that for the strong shock nilri\ ~ 4 for Y = 5/3 and nijn\ <~ 7 for 
r = 4/3. In the Schwarzschild black hole case, for example, we obtain a similar tendency for the maximum 
values; that is, n^jrix ~ 4 for Y — 5/3 and n^lrix > 8 for Y = 1.4 (see Figs. 3c and 4b). The compression 
ratio also has a maximum £ max w 3.6 for a = and £ max > 4 for a = 0.95m (Fig. 3a). Compared with the 
n 2 /ni vs v\ diagram, £(v[) is almost constant for the strong shocks (i.e., a > 0). This is due to the effect of 
factor u\/u\ (see, cq. (56)). From Figures 3a and 4b, it may appear that a strong shock is possible for very 
low accretion speed (for a > —0.2m). Note that in our calculation, we set r\ = constant (77 oc nu r ^ 0), so 
means rii > 1. When d[<1, the postshock flow must also satisfy *C 1, and then n 2 3> 1. Thus, 
as a result of numerical calculation, we obtain a finite value for 77.2/ni in the limit of v[ — > 0. 

Figures 5a and 5b show the £ vs. vf and vf vs. v{ relations, where v$ = (— g^/ p w )(^l — w) is the toroidal 
3-velocity seen by a ZAMO. Figure 5b shows that when its radial velocity is very low, such a flow does possess 
sufficiently large toroidal velocity (especially the a = 0.95m case). The negative vf is due to the effect of 
dragging on the plasma by the counter-rotating black hole. The toroidal velocity vf is explicitly related to 
the toroidal magnetic field (B$)i, and therefore it can be also related to the magnetization parameter a\ 
(see Fig. 4c). The difference of o\ between the corotating black hole case and the counter- rotating black hole 
case is caused by this spin-dependence of the toroidal motion. That is, due to the dragging effects, \B^\\ is 
forced to decrease for a — 0.95m, while it is forced to increase by stretching toward the (—(^-direction for 

the a = —0.4m case. We also obtain the result that (vf — vf) — > and the magnetic field line is tangent to 
the shock surface when v^ —* 0. Figure 5c shows the £ vs. r s h relation. For a magnetosphere with a large 
spin the shock front is located near the event horizon, and a strong shock is generated. Along each curve 
the shock strength increases slightly when the shock location shifts inward from the vicinity of the plasma 
source, but near the Alfven point it decreases considerably. 

We extended the work by AC88 on the MHD shock conditions for special relativistic jets, by deriving 
corresponding equations for general relativistic accretion flows in the magnetosphere of a rotating black hole. 
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Our results refer to the slow magnetosonic shock for accretion flows very close to a black hole where the 
general relativistic effects are strong, while the AC88 results refer to the fast magnetosonic shock for outflows 
at long distances away from the hole. The one apparent difference between our results and those by AC88 
is that their r (our £) and II — > sharply as their (3 — ► 0, while our £ and II do not drop off sharply as v\ 
— ► (see our Figures 3a and 3b). (In the rest of the velocity range, we find that the general trend of our 
results and theirs is quite similar, especially for positive a.) Another difference is the velocity dependence 
of q. For the fast magnetosonic shock in the case of AC88, there is a lower-limit threshold for the upstream 
velocity v\ to satisfy the requirement of q > 1. When this lower-limit velocity is approached (where q = 1), 
the compression ratio decreases, and below this velocity there is no fast magnetosonic shock solutions. On 
the other hand, for the slow magnetosonic shock in our case, the condition q < 1 is required, and moreover 
there is not such a threshold for the lower limit to the upstream velocity and £ does not drop off. For v\ — > 
we obtain q = 1, and for the switch-off shock (if present) we obtain q = (see Fig. 4a). The switch-off shock 
gives the upper limit to the upstream flow velocity, which is just the Alfven wave velocity. Note that in the 
rotating black hole case the velocity range that allows to the slow magnetosonic shock is restricted by some 
other factors — for instance, due to the dragging effect. 

The general behavior of T and 9 with velocity for all spin parameters a as exhibited in Figures 3c 
and 3d are essentially the same as Figure 6 of AC88 for hydrodynamic shocks with no plasma rotation. As 
mentioned by AC88, we also find that the polytropic index for the postshock plasma stays closer to 5/3 than 
to 4/3, except for extremely high Lorentz factors; here we should note that, for small v\, T remains ^5/3 
even when vf dominates the flow velocity. 

We note that our equations derived in Section 3.2 reduce to the corresponding equations in AC88 in the 
limit of no plasma rotation and weak gravity (Qi = 0, m = and a = 0). With more direct comparison with 
AC88 in mind, we present in Figure 6 the result for the asymptotic case of a = 0, Oi = 0, slowly rotating 
magnetosphere Slf « 1, and large distance, closer to the situation investigated by AC88. Figure 6a shows 
the compression ratio vs. v\ relation, while in Figure 6b the adiabatic exponent and temperature parameter 
are shown as a function of v\. The chosen conserved quantities are given in Table 3. These parameters 
appear to put us in the small b' and small 8' regime (small a) of AC88. Figures 6a and 6b show that 
although the parameters chosen here and AC88 are somewhat different, these curves are essentially identical 
to figures 1 and 6 of AC88. (Note that figure 6 in AC88 shows the hydrodynamical shock case.) 

Though our results of slow magnetosonic shocks for a > are very similar to the AC88 results obtained 
for the fast magnetosonic shocks, the local physical properties are obviously different between the fast and 
slow magnetosonic shocks. Here, we should note that our results include the rotational effects of magnetic 
field lines, which arc leading the black hole toward the rotational direction. So this apparent similarity is 
caused essentially by the rotation of the magnetic field lines. When < Oi < Of, we find that Oi < 2 < Of- 
This means that in our slow magnetosonic shock case the trajectory of a fluid element seen by a distant 
observer becomes very similar to the fast magnetosonic shock case without rotation, as shown in Figure 7a. 
On the other hand, in the shock frame rotating with Of, the trajectory of a fluid element in our case is 
coincident with the magnetic field pattern (see Fig. 7b). So, some shock properties of the slow magnetosonic 
shock with Oi > behave like the fast magnetosonic shock without rotation (especially the hydrodynamical 
behavior). This is not essentially due to the general relativistic effects, but the existence of the slowly 
rotating black hole (0 < < Of), or a counter-rotating black hole (fipWH < 0), is responsible for such a 
situation (with the magnetic field lines leading the black hole in rotation) in the rotating magnetosphere. 
When Oi < 0, which state is easily realized for a rapidly counter-rotating black hole, the downstream plasma 
angular velocity 2 is in the range of Oi < 2 < or < 2 < Of, and the trajectory of the fluid is quite 
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different from the fast magnetosonic shock case, and then the shock behavior is also different. Also, for a 
rapidly rotating black hole with < Qp < and the magnetic field lines trailing the black hole in the 
rotational direction, this similarity would be lost except for the very strong shock cases. Note that the 
dominant energy dissipation and accompanying radiation processes in the postshock hot plasma would be 
different between the corotating slow magnetosonic shock and the fast magnetosonic shock without rotation, 
even if some of their properties are very similar. 

Finally in Figure 8 we show our preliminary results for the upstream (solid curve) and downstream 
(dashed curve) magnetization parameters, o\ and (72, plotted as a function of the shock location, r s h, for 
a representative spin of a = 0.7m. The conserved quantities chosen for this spin parameter are given in 
Table 2. From the energy momentum conservation at the shock front with equation (51), the postshock 
magnetization parameter <J2 is expressed as 02 = C~ 1(J i + (C _1 — !)• The decreasing a across the shock 
means that £ > 1. Although, at the shock front, kinetic energy of the preshock flow obviously converts to 
thermal energy of the postshock flow, the electromagnetic energy also converts to thermal energy at the same 
shock front. At the location of r s h — > r\ (when the switch-off shock takes place), the energy conversion is 
maximum. The detailed transitions of each of the energy components and the discussion of the implications 
will be presented in Paper II (Rilett et al., in preparation). 

5. Summary and Concluding Remarks 

The main purpose of this paper is to explore the general relativistic version of MHD shock conditions 
for accreting plasma in a magnetosphere of a rotating black hole. Here, we introduce the magnetization pa- 
rameter seen by ZAMO, because the shock conditions are related to the local plasma and the magnetosphere 
is dragged by the rotating black hole. Then, wc can obtain similar expressions as for the special relativistic 
case, though the formulae include the Kerr metric. We have seen that the compression ratio £ is the solution 
of a polynomial of eighth degree. 

Although the jump condition is of course determined by the local physics, for a distant observer, the 
shock properties depend also on the gravitational field structure (i.e., the gravitational redshift factor and 
the dragging effect, etc); furthermore, the preshock plasma flow parameters are also affected by the gravity 
and its spin. So, when we observe some physical properties of the magnetosonic shock, which may be related 
to some efficient radiation process, we should understand the gravitational field structure in addition to the 
properties of magnetized plasma flow. 

In this paper our sample solutions are presented for the special case of slow magnetosonic shocks in 
plasmas accreting in the equatorial plane with the radial magnetic field configuration. For slow magne- 
tosonic shocks our results are compared with those by AC88. It is confirmed that, though AC88 discussed 
fast magnetosonic shocks, the general agreement is excellent in the limit of no plasma rotation and weak 
gravity. However, with the introduction of additional factors, gravity and rotation, we find that the situation 
becomes far more complicated. For example, in a black hole magnetosphere, two Alfven points and two slow 
magnetosonic points can be located inside the separation point r = r sp , where the gravitational force and the 
centrifugal force are balanced for plasma with zero poloidal velocity and the radius is defined by (a') sp = 0; 
the prime means the derivative along a stream line. In this paper we have treated only cold MHD accreting 
flow passing through the inner Alfven point. However, the MHD accretion flow passing through the outer 
Alfven point is also possible. In both cases the injection point can be located near (inside) the separation 
point. When we consider shock formation where the postshock flow passes through the outer Alfven point, 
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we would find such a flow as an additional branch in the figures such as those presented in this paper. 
Such a new branch would appear outside the previous branch, and it of course must be located between the 
outer Alfven point and the injection point. However, the postshock accretion flow may pass through the 
inner Alfven point. Furthermore, to obtain the accretion flow passing through the fast magnetosonic point 
twice by the shock formation, we must more carefully find the acceptable parameter ranges of conserved 
quantities. This will be our next step. In non-relativistic flow solutions, Chakrabarti (1990b) discussed the 
slow and fast magnetosonic point conditions and shock conditions in winds, where multi-magnetosonic wind 
solutions were demonstrated. In the future, we will extend our results by treating a multi-magnetosonic flow 
conditions with shock formation in the Kerr geometry. 

As a preliminary example of how the transition of various parameters behave across the shock front, 
in Figure 8 we showed the prcshock and postshock magnetization parameters. The detailed results for the 
transition behavior of various energy components and other parameters will be presented in Paper II (Rilctt 
ct al., in preparation). Here we only comment that both fluid kinetic energy and electromagnetic energy 
convert to thermal energy at the shock location. 

When we consider accretion through high-latitude regions above the event horizon near the hole's 
rotational axis (see Nitta, Takahashi & Tomimatsu (1991); Tomimatsu & Takahashi (2001)), we can expect 
as an X-ray source hot ringed area around the axis. Our extended future work will include the application of 
our current work to realistic models of a black hole magnctosphere as the central engine of accretion-powered 
AGN. 

We are grateful to Akira Tomimatsu and Hideyo Kunieda for useful discussions. S.T. thanks the Yamada 
Foundation for financial support through the visiting program. 
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r/m 

Fig. 1. — A preshock cold trans- fast MHD accretion solution (rj = r/p; bold curve with arrow). The solution 
with zero- velocity at the event horizon r = rjj is unphysical. The other solutions (77 ^ tjf; thin curves) are 
also unphysical. The radii of r = rp , r = tl , r = ta and r = n n j are the locations of the fast magnetosonic 
point, the light surface, the Alfven point and the injection point, respectively. A slow magnetosonic shock 
would be possible somewhere between the plasma injection radius and the Alfven radius. 




Fig. 2. — Accretion flows with shock front. It is assumed that in the poloidal plane the downstream flow is 
radial and the shock front distributes perpendicular to the the downstream flow direction. In a general case, 
the magnetic field has a non-zero toroidal component (marked by <g>) because of the plasma rotation. 
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Fig. 3. — (a) The shock frame compression ratio £, (b) thermal pressure II, (c) polytropic index T, and 
(d) temperature parameter O, respectively, as a function of v\, the radial velocity of the infalling plasma, 
for variable spin parameters (a/m = 0.95,0.3,0.0,-0.2,-0.4). The parameters II, T and O refer to the 
downstream values. The shocks at the switch-off points are shown by small spots. The sets of field aligned 
parameters are given in Table 2. 



- 18 - 



a) 



1.0 



. ! 



. 6 



0.4 



0.2 



0.0 





^ -0.4 












\\ 




0.95\ a /m=0.0 \\ 




\ l| 
\ 1 




\ " 



0.0 0.2 0.4 0.6 

1^ 



l . o 



(b; 



10 . 



n 1 



8 . 



6.0 



2 . 





t 

a/m=0.3 1 






1 /» 

//o.o 




0.95/ y // /'\ 




/-OA 

/ 







0.0 0.2 0.4 0.( 



1.0 



(c) 



o . f 



. 6 



0.4 



0.2 



0.0 



/- 0.4 

i 

i 

! / 

/ y-o.i 

a/m=0.0 



0.3 



0.95 




OA 



1 . 



Fig. 4. — Various parameters are shown as a function of v\. (a) The ratio of the downstream toroidal 
magnetic field to upstream toroidal magnetic field q, (b) the ratio of the downstream number density to 
upstream number density, and (c) the upstream magnetization parameter o\. The chosen parameter sets 
are the same as in Fig. 3. 



- 19 - 



(a) 



4.0 



3.5 



3.0 



2.5 



2.0 



OA 



0.4 - 





! J a/m=0.0 

HI 


0.95 


-0.4 


i 1 \ 

1 1 N 




T\ 


1 i -0.2 




. ! 







-0.2 - 



-0.2 0.0 


0.2 0.4 


. 0. 






(b) 






°- 95 






-_0.3 






0.0 o 












-0.2 




a /w- O. I 





0.0 0.2 0.4 0.6 0.8 1.0 



V 



(c) 



4.0 



3.5 



3.0 



2.5 - 



2.0 



0.95 



/ 

/ /' 

V 

11 



0.3 



a/m=0.0 



-0.2 



\ -0.4 

\ 

\ 

\ 

3 4 5 

r sh /m 



Fig. 5. — (a) The compression ratio as a function of vf, (b) the toroidal velocity vf as a function of v{, and 
(c) the compression ratio £ as a function of the shock location r s h . The parameter sets are the same as in 
Fig. 3. 
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Fig. 6. — (a) The shock frame compression ratio £, and (b) polytropic index T and temperature parameter 
0, respectively, as a function of the radial preshock velocity v\ in the limit of no plasma rotation and weak 
gravity (f2i = 0, m = 0, a = 0). The set of field aligned parameters is given in Table 3. 



- 21 - 




Fig. 7. — The slow magnetosonic shock seen by (a) a distant observer in the (r, </>)-plane and (b) a shock 
frame observer in the (r, 0)-planc, which is corotating with the magnetic field line. For a slowly rotating 
black hole case (wh < Of), the magnetic field lines are leading the black hole toward the rotational direction, 
and the fluid trajectory across the shock front looks like the fast magnetosonic shock case; and then some 
physical properties of the slow magnetosonic shock is very similar to that of the fast magnetosonic shock. 
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Fig. 8. — The upstream (solid curve) and downstream (dashed curve) magnetization parameters, u\ and 02, 
as a function of the shock location, for a = 0.7m (see Table 2 for chosen flow parameters). 
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Conserved 
Parameter 


Physical description 


Definition 


n F 


the angular velocity of the magnetic field lines 


Eq.(8) 


E 


the total energy seen by a distant observer 


Eq.(17) 


L 


the total angular momentum seen by a distant observer 


Eq.(18) 


V 


the particle number flux per magnetic flux tube 


Eq.(19) 



Table 1: The table of conserved flow parameters. The mathematical definitions are expressed in the text 
(see, Bekenstein & Oron (1978); Camenzind (1986a)). 



a/m 




rnQp 


rn/m 


fij/m 


r^/m 


r inj /m 


r sp /m 


0.95 


2.470 


0.3691 


1.312 


1.314 


1.335 


1.439 


1.457 


0.70 


3.137 


0.2214 


1.714 


1.718 


1.757 


2.337 


2.442 


0.30 


3.950 


0.1096 


1.954 


1.964 


2.022 


3.694 


4.270 


0.00 


5.073 


0.0516 


2.000 


2.022 


2.093 


5.508 


7.256 


-0.20 


7.332 


0.0278 


1.980 


2.029 


2.102 


6.197 


8.504 


-0.40 


10.233 


0.0062 


1.916 


2.010 


2.090 


21.406 


29.609 



Table 2: The parameter sets for cold trans-fast MHD accretion solutions, where Ej [i\ — 1.006, \i\r\ = 0.041 
and fii — m p (for a cold flow). 



a/m 


E/fti 








r L /m 


r A /m 


n-n.il m 


0.0 


1.00024 


15.000 


0.0019 


0.0005 


2.00003 


2.001 


44.6 



Table 3: The parameter sets for cold trans- fast MHD accretion solutions. 



